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Abstract 

Neutral Higgs-boson pair production provides the possibility of studying the 
trilinear Higgs couplings at future high-energy colliders. We present the QCD 
corrections to the gluon-initiated processes in the limit of a heavy top quark 
in the loops and the Drell-Yan-like pair production of scalar and pseudoscalar 
Higgs particles. The pp cross sections are discussed for LHC energies within 
the Standard Model and its minimal supersymmetric extension. The QCD 
corrections are large, enhancing the total cross sections significantly. 
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1 Introduction 

The Higgs mechanism is a cornerstone of the Standard Model (SM) and its supersymmet- 
ric extensions. Thus, the search for Higgs bosons is one of the most important endeavours 
at future high-energy experiments. In the SM one Higgs doublet, has to be introduced 
in order to break the electroweak symmetry, leading to the existence of one elementary 
Higgs boson, H |IJ . The scalar sector of the SM is uniquely fixed by the vacuum expecta- 
tion value v of the Higgs doublet and the mass M# of the physical Higgs boson . Once a 
Higgs particle is found, it is necessary to investigate its properties in order to reconstruct 
the Higgs potential and to verify that it is indeed the SM Higgs boson. A first step in this 
direction is the measurement of the trilinear self-couplings, which are uniquely specified 
by the scalar potential 



The parameter A defines the strength of the Higgs self-interactions. In the SM it is given 
by A = Mfj/(2v 2 ). At tree level, A can only be probed through multiple Higgs-boson 
interactions and there are, at present, no direct experimental limits on A. In extensions 
of the SM, such as models with an extended scalar sector, with composite particles or 
with supersymmetric partners, the self-couplings of the Higgs boson may be significantly 
different from the SM predictions. The limits that may be obtained on the trilinear self- 
coupling of the Higgs boson at the LHC and the impact of QCD corrections represent a 
particular topic of this paper. 

Since the minimal supersymmetric extension of the Standard Model (MSSM) requires 
the introduction of two Higgs doublets in order to preserve supersymmetry, there are five 
elementary Higgs particles, two CP-even (h, one CP-odd (A) and two charged ones 
(H^). This leads to a large variety of self-interactions among them. At lowest order all 
couplings and masses of the MSSM Higgs sector are fixed by two independent input param- 
eters, which are generally chosen as tg/3 = V2/V1, the ratio of the two vacuum expectation 
values vi t 2, and the pseudoscalar Higgs-boson mass Ma- The self-interactions among the 
Higgs bosons (at lowest order) are given in terms of the electroweak gauge couplings, 
tg/3 and Ma, and may be quite different from the ones of the SM, which are governed 
by the parameter A. Higher-order corrections to the MSSM Higgs sector turn out to be 
important owing to the large top-quark mass rrtt |3|-[5|]. They increase the upper bound on 
the light scalar Higgs mass M^ from the Z-boson mass Mz to about 130 GeV, along with 
altering the Higgs-boson self-couplings with contributions proportional to G pm\ j M\. 

Higgs-boson pairs can be produced by several mechanisms at hadron colliders: 

• Higgs-strahlung W*/Z* -> fafcW/Z [§, 

• vector-boson fusion WW, ZZ — > (fiifo 0, 

• associated production Z* — > Ah, AH, 

1 We have taken care that no confusion can arise from using the same symbol H for the SM and the 
heavy CP-even MSSM Higgs particle. 
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Higgs radiation off top and bottom quarks gg, qq — > QQ4>\(j)2 
gluon-gluon collisions gg — > 0102 R |10 ■ 



At the LHC, gluon fusion is the dominant source of Higgs-boson pairs, although in some 



regions of the MSSM parameter space, vector-boson fusion |TT| can be important. Note, 
however, that gg — > HA represents an exceptional case, since this channel is suppressed 
with respect to the Drell-Yan-like process qq —>■ Z* — > HA, so that it will be very 
difficult to separate the gluon-fusion process in this case experimentally. The gluon- 
fusion process gg — > hA, on the other hand, is competitive with the Drell-Yan-like process 
qq -> Z* -> /iA 

In this paper we present the QCD corrections to the Drell-Yan-like production and 
the gluon-gluon collision processes. The gluon-fusion processes are, in the SM, mediated 
by triangle and box loops of top and bottom quarks; in the SM, the contributions of the 
bottom quark can always be neglected. In the MSSM, the squark contributions will be 
suppressed if the squarks are heavier than ~ 400 GeV, and, for small tg/3, the top-quark 
loops dominate the gg cross sections. The QCD corrections to the gluon-fusion processes 
have been obtained in the limit of a heavy top quark by means of low-energy theorems 
and also by explicitly expanding all relevant one- and two-loop diagrams. The results 
are expected to be valid for small tg/3 in the MSSM and below the tt threshold of the 
top-quark loops in both the MSSM and the SM, since in this regime effects of a finite 
top-quark mass are expected to be small. In the case of single-Higgs production the same 
procedure reproduces the known exact result for the NLO cross section within 5%, for 
Higgs-boson masses below 2m t [12|. The considered QCD corrections are important in 



the process of extracting limits on the Higgs-boson self-couplings reliably. 

The paper is organized as follows. In the next section the low-energy theorems for 
the interactions of gluons with light Higgs bosons will be reviewed, and the relevant 
interactions in the heavy-quark limit will be constructed. In Section |^ the details of the 
calculation will be described, and in Section [| we present the results for the SM and 
MSSM Higgs bosons. In Section |5| we give some conclusions. 



2 Low-Energy Theorems 

In the low-energy limit of vanishing Higgs four-momentum, the Higgs-field operator acts 
as a constant field. In this limit it is possible to derive an effective Lagrangian for the 
interactions of the Higgs bosons with gauge bosons, which is valid for light Higgs bosons. 
This effective Lagrangian has been successfully used to compute the QCD corrections to 
a number of processes, in particular to single-Higgs production from gluon fusion at the 



LHC [[13|, |i4| . In this case, the result of using the low-energy theorems has been shown 
to agree with the exact two-loop calculation to better than 10% even for M H as large 
as 1 TeV. This lends legitimacy to our use of the low-energy theorems to compute QCD 
corrections to multiple Higgs-boson production via gluon fusion. 
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In the limit of vanishing Higgs four-momentum the entire interaction of the scalar 
Higgs particles if, with a heavy quark Q can be generated by the substitution 



771} 



<+V° Q Y,9h H i 



(2) 



in the Lagrangian of a heavy quark of bare mass rrin ]Tj|, where t]q = ttiq/v denotes the 
bare SM Yukawa coupling, which must not be included in the substitution. The symbol 
g l q is the relative strength of the heavy-quark Yukawa coupling, 



£ Yu k = ~Vq 9qQQHi. 



(3) 



In the SM, we have g l q = 1. The expressions for g l q in the MSSM are given in Ref. [|14 



At higher orders this substitution has to be performed for the unrenormalized parameters 
In the following we restrict our analysis to the top-quark contributions. At 
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next-to-leading order (NLO) the effective interaction between several scalar Higgs fields 
and gluons can be obtained from the radiatively corrected effective Lagrangian of gluon 
fields, 

1 



99 



1 + 14(0) 



(4) 



with n* g (0) denoting the top-quark contribution to the unrenormalized gluon vacuum 
polarization at zero-momentum transfer. At two-loop order, we have 
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(5) 

where the strong coupling constant includes five light flavours. This means that the 
top-quark contribution to the running of a s has been subtracted at vanishing momentum 
transfer. Hereafter, we drop the superscript 5 on a s . Performing the substitution Eq. (fj) 
and renormalizing the bare top mass m" via 



m. 



i-^r(i + e) {^4- 



where m t denotes the pole mass, we end up with the NLO result 
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The interaction of even numbers of pseudoscalar Higgs bosons with gluons can be 
determined analogously from Eq. (|) by substituting [16|] 



K) 2 - K) 2 + (gfvtA) 2 , 



leading to 
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The case of odd numbers of pseudoscalar Higgs bosons can be derived from the ABJ 



anomaly |L7j in the divergence of the axial vector current |14| , [lq , pL8| . The interactions 
that are relevant in our case are0 

r — n A as r ,a i iv n a ^ 

L-Agg ~ Q t ~ , 

£ah 99 = ~9tg?^G a ^ v G a ^—^-, (10) 
where G a = ^e^vapG 0,013 denotes the dual gluon field-strength tensor. Owing to the 



Adler-Bardeen theorem jnj there are no higher-order corrections to the effective La- 
grangians involving odd numbers of pseudoscalar Higgs bosons. 

Figure [l| summarizes the Feynman rules for the effective interactions between two 
gluons and one or two Higgs bosons; the rules can be read off from Eqs. (|7|)- (|10D . These 
Feynman rules can now be used to compute Higgs interactions beyond the lowest order. 
We recall that there is no contribution of light quarks (which are considered to be massless) 
to the effective couplings, but note that light-quark loops have to be included when the 
Higgs bosons do not directly couple to the quark loops. Such contributions arise, in 
particular, in gg — > Z* — > hA, HA and cannot be obtained from the low-energy theorems. 



3 QCD Corrections 

3.1 Gluon fusion: basic definitions 

At leading order (LO) neutral-Higgs pair production via gluon fusion is mediated by 
triangle and box diagrams of heavy quarks, as exemplified in Fig. |2|. In the heavy-quark 
limit, the fermion triangles and boxes can be replaced by the effective vertices of Fig. [I]. 
Throughout this analysis, we choose the squark masses to be 1 TeV so that squark-loop 
contributions can be neglected in the MSSM case. Generically the partonic LO cross 



section can be expressed as JTU 



>W2 



dt 



^3 I 
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256(2tt) 

The Mandelstam variables for the parton process are given by 
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m 2 ~ v HQ 2 i m \i m i) cos $ 



m i ~ m l + \J ^(Q 2 i m ii m i) cos 



(12) 



2 Note that in the earlier Refs. JTJ, a factor of 1/2 is missing in the effective Lagrangians for 

the single pseudoscalar Higgs-boson coupling to gluons. 
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-ita^g^g? 2 {-<T(*i ■ h) + k»k%} (1 + ^) 




i5 a^{gtf {-9T(kx ■ h) + %k$} (1 + ^) 




Figure 1: Feynman rules for the effective interactions of Higgs bosons with gluons in the 
heavy-quark limit, including NLO corrections. 
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Figure 2: Generic diagrams describing neutral Higgs-boson pair production in gluon- 
gluon collisions (<p,<pi = h,H,A). 



where 6 is the scattering angle in the partonic cm. system with invariant mass Q, and 

z) = (x - y - z) 2 - Ayz. (13) 

The integration limits 

i± = - 1 \q 2 -ml-mlT J HQ 2 , ml m|)l (14) 



in Eq. (JlTJ) correspond to cos# = ±1. The scale parameter /x is the renormalization scale. 
The complete dependence on the fermion masses is contained in the functions Fa, F d , and 
G n - The full expressions of the form factors Fa, F n , G n , including the exact dependence 
on the fermion masses, can be found in Ref. [IU|. 

The couplings Ca and C a and the form factors Fa, F n ,G n in the heavy-quark limit 
are given by: 



(i) SMi 
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with the trilinear coupling Xhhh = 3M^/Mf. 
;ii) MSSM: 

4>i4>2 are generically defined as 



The couplings for the processes gg 
h, H, A) 



ci 



A 



M 2 



4>i4>24> 



s-M 2 + iM 6 r 



-9t 



4>i „4>2 



Cn = gTg\ 



(16) 



where denotes the Higgs particles of the s-channel contributions. The trilinear 
couplings Xfafatj) and the normalized Yukawa couplings gf can be found in Ref. . 
The individual expressions in the heavy-quark limit can be summarized as: 
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0i0 2 = hh, hH, HH: 
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(17) 



W2 - 



C A - C£ + C& Cl - X ZAh , ZAH J_~M^M^ at ' 

Ft - 1, it-^^(M/U-Mi), F D --1, 

z 

G n -> 0, 

where a 4 = 1 denotes the axial charge of the top quark. 
A4: 



C A 


= cl 


+ C7f, 


F a 
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3' 



(19) 



It should be noted that owing to the Ward identities for the Zgg vertex only the pseu- 
doscalar Goldstone component of the Z bosons contributes to in the case of HA and 
HA production in the heavy-quark limit. The anomaly contributions of the top and 
bottom quarks cancel. 

The QCD corrections consist of two-loop virtual corrections, generated by gluon ex- 
change between the quark lines and/or external gluons, and one- loop real corrections with 
an additional gluon or light quark in the final state. We have evaluated the QCD correc- 
tions in the heavy-quark limit by means of two different methods: (i) using the effective 
couplings based on the low-energy theorems, as presented in the previous section, and (ii) 
explicitly expanding all relevant one- and two-loop diagrams in the inverse heavy-quark 
mass. In the following we shall describe the details of both approaches. 



3.2 Low-energy theorems 

We are now in a position to compute the NLO corrections to Higgs-boson pair production. 
Typical effective diagrams contributing to the virtual and real corrections are presented in 
Fig. |3[ Adopting the Feynman rules of Fig. [I] for the effective interactions, the calculation 
has been carried out in dimensional regularization with n = 4 — 2e dimensions. The strong 
coupling has been renormalized in the MS scheme including five light-quark flavours, i.e. 
decoupling the top quark in the running of a s . After summing the virtual and real 
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Figure 3: Typical effective diagrams contributing to the (a) virtual and (b) real corrections 
to neutral Higgs-boson pair production. 



corrections the infrared singularities cancel. However, collinear initial-state singularities 
are left over in the partonic cross sections. Those divergences have been absorbed into 
the NLO parton densities, defined in the MS scheme with five light-quark flavours. We 
end up with finite results, which can be cast into the form 



0-KSLO{PP -> 0102 + X) = (X L o + A(T v i r t + A(T gg + Aa gq + Aa qq , 

with the individual contributions 
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The coefficient C for the virtual corrections reads 
2 , 33-2N F l /i 2 



// + dt \c 2 C n (C A F A + CaF u ) + c 3 §|(Q 2 - m\ - m 2 )C 2 G D 

+ 3fJe — ^ -, >- (22) 

J + dt{\C A F A + C n F n | 2 + |C D G n | 2 } 



where 



_ (m 1 +m 2 ) 2 2 _ {t - m\){u - m\) , 

r — > Pt — Q2 m l- ^ 6 > 



The objects P gg (z) , P gq (z) denote the Altarelli-Parisi splitting functions [^0 

PM " 1 

pm = ^ 1 v ; - , (24) 

where Np = 5 in our case. The factorization scale of the parton-parton luminosities 
dD 3 jdr is denoted by M. The coefficients q for the individual final-state Higgs bosons 
0102 are given by 




II 4 4 

9' ° 3 ~ ~9 ; 



'102 = hh,hH,HH : ci = — , c 2 = -, c 3 = --, 



L ^ rr^ a 2 2 t~U 

h<p2 — hA,HA: ci = 6, c 2 = -, c 3 



3' ~° 3 g 2 -m 2 -m|' 
0102 = AA : Ci = y, c 2 = -1, c 3 = -1. 

(25) 

In order to improve the validity of our results, we have inserted the full expressions for the 
form factors F A , F n and G a in Eqs. ([□]) and (|22]) , i.e. including the exact dependence on 
the fermion masses. This procedure is reasonable since the QCD corrections are dominated 
by soft and collinear gluon radiation, which do not resolve the Higgs couplings to gluons, 



analogously to single-Higgs production via gluon fusion |L2 . 

A few remarks on the s-channel Z-boson exchange in hA, HA production are in order. 
For the virtual corrections, the factorization of the NLO corrections into the LO form 
factors and a universal correction factor is exact for m t — > oo and = 0. This is due to 
the fact that only the pseudoscalar Goldstone component of the Z boson contributes as 
in LO, i.e. the QCD corrections coincide with the one to s-channel pseudoscalar Higgs- 
boson exchange. For the real corrections, the factorization is not exact, but the applied 
correction factor correctly includes the dominant contributions, which are caused by soft 
and collinear gluon radiation. Additional infrared- and collinear-finite contributions, e.g. 
originating from Zggg box corrections in gg — > g + Z* g + hA, HA processes, are 
expected to be small, since they do not exhibit large contributions from soft and collinear 
gluon radiation. They are neglected in our analysis. 
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3.3 Explicit expansion of the gluon-fusion diagrams 

We have derived the above results also by explicitly performing the heavy-mass expansion 
of the corresponding one- and two-loop Feynman diagrams. The amplitudes for the indi- 
vidual diagrams have been automatically generated using the package FeynArts |21] . The 
asymptotic expansion of the individual amplitudes in the heavy top-quark mass is carried 
out directly in the integrand, i.e. before the integration over the momentum space. We 



employ the general algorithm of Ref. |22[ (see also Ref. |23] and references therein) for the 
asymptotic expansion of Feynman diagrams in dimensional regularization. This method 
expresses the coefficients of the expansion in terms of simpler diagrams. At the one-loop 
level, this procedure leads to simple one-loop vacuum integrals only. At the two-loop level, 
we get two-loop vacuum integrals and products of one-loop vacuum integrals and massless 
one-loop integrals with non- vanishing external momenta. The analytical calculation of all 
those integrals is straightforward when using the Feynman-parameter technique. Since 
the employed strategy leads to a very large number of terms in intermediate steps, and 
since each step is algorithmic, we have fully automatized the calculation in Mathematica 



24}| . In the following we sketch the single steps of the calculation and give the results for 
the basic integrals. 

The general algorithm for the asymptotic expansion of any given Feynman graph F 
in the limit Mi — *> oo for some internal masses Mi can easily be summarized. Denoting 
the corresponding Feynman amplitude by F? and the corresponding integrand by It, the 
large-mass expansion reads 

= J(Ud n n)lr ^ E J (l[d n qi) Ir h r plmi I,, (26) 

where qi are the integration momenta. The sum on the r.h.s. runs over all subgraphs 7 
of T which contain all propagators with the heavy masses Mj and which are irreducible 
with respect to those lines of 7 that carry light masses mi. The integrand of the subgraph 
7 is denoted by J 7 . The reduced graph r/7 results from V upon shrinking 7 to a point, 
and the integrand Jr/ 7 is defined such that Jr = / 7 /r/ 7 - The symbol T p i mi represents an 
operator that replaces the integrand J 7 by its Taylor series in the expansion parameters 
pj and m>i, where pj are the external momenta of the subgraph 7. Therefore, Eq. (^) 
expresses the original integral Fp in terms of an infinite sum over simpler integrals. For any 
given power £ a , this sum contains only a finite number of terms that are non-vanishing 
in (-Fr|ju i -+$M i )/C a after the scaling limit £ — ► 00 is taken. These terms can easily be 
determined by power counting. This general strategy for the heavy-mass expansion will 
become more transparent when we inspect in more detail the different types of graphs 
that are relevant in our case. 

We start by considering the relevant one-loop integrals in the limit m t — > 00. They 
contain only top-quark propagators in the loop, both for the LO calculation and for 
the real NLO corrections. According to the algorithm (|26f) , there is only one relevant 
subgraph 7, namely the subgraph 7i oop containing only the propagators of lines inside the 
loop. If T is irreducible, we have T = 7i oop - The Taylor-expansion operator T replaces 
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each propagator P(q — p,m t ) by 



P(q - p, m t ) = [(q - p) 2 - m 2 ]' 1 = ]T(g 2 - m 2 )- 1 ~ l (2qp - p 2 )\ (27) 

1=0 

where q denotes the integration momentum, and p is any combination of external mo- 
menta. These replacements express each one-loop diagram by a sum of terms containing 
one-loop vacuum integrals 

«,.(»^.) = ^/A^. (28) 

The terms that are non-vanishing in the heavy-mass limit can be determined by simple 
power counting, since an explicit factor of m t and the integration momentum q contribute 
to the scaling factor in m t exactly in the same way. All non-vanishing vacuum tensor 
integrals Q28D can be decomposed into terms that are products of metric tensors g^ and 
coefficient factors. The coefficients for the different covariants, which span the whole 
tensor, can be algebraically expressed in terms of scalar vacuum integrals. This algebraic 
reduction, which proceeds recursively in the tensor rank, is standard. The trick is to 
contract the equation that expresses the integral in terms of covariants with a set of some 
suitably chosen covariants. On the side of the integral, this leads to integrals that are 
already known; on the other side of this equation, one gets linear combinations for the 
tensor coefficients. The coefficients are obtained by inverting a system of such linear 
equations for the coefficients. For the case of one-loop tensor vacuum integrals, at most a 
single covariant structure contributes, namely the totally symmetric tensor built of metric 
tensors g^. Only scalar one-loop integrals, i.e. the ones of (|28| ) with R = 0, have to be 
computed explicitly. A simple calculation yields 

V^im-mx) = (-l)"HV)^m"- 2ni ^tr^i. (29) 

r(n x ) 

At the two-loop level, there are two basically different types of diagrams. The first 
type contains two independent top-quark loops. Such diagrams do not lead to genuine 
two-loop integrals and can be treated like the one-loop diagrams above. Topologically 
those diagrams are represented by the third graph of Fig. [| The second type of graphs 
is formed by the genuine two-loop diagrams. Each of those diagrams contains a closed 
top-quark loop and one, two, or three internal gluon lines. Typical box diagrams are 
shown in Fig. f|. For all such genuine two-loop diagrams there are two subgraphs 7 that 
are relevant in the expansion ( ffif ) for m t —>■ 00. The first subgraph is the diagram 7i oop 
built of all lines inside the loops, the second is given by the closed top-quark loop. 

First we consider the case 7 = 7i oop - The Taylor expansion of the integrand J 7loo 
involves the consistent expansion of each propagator about the external momenta of the 
process. This means that each propagator P(q — p, m) is replaced by 

00 

P{q-p,m) = [(q-p) 2 - m 2 }' 1 = ^(g 2 - m 2 )~" 1-, (2gp -p 2 )', m = 0,m t , (30) 

1=0 
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Figure 4: Typical two-loop box diagrams with one, two, or three internal gluons. 



where q is a linear combination of the two integration momenta qi and q 2 , and p consists of 
external momenta pi. This replacement leads to two- loop vacuum integrals of the general 
form 



v t£ ) ... m -v 1 ...v s { n ii n 2, n 3 ; mi, m 2 , m 3 ) 



(27T/X) 



27T^ J t7T z J 



qi 



■ qi,n R 92,1/1 • • • 92. 



{ql-ml)^{ql 



92 J 



(31) 



with the actual mass insertions m\ = and m 2 = m 3 — m t . In the power counting, which 
determines the non- vanishing terms for m t — > oo, factors of the integration momenta 
5i and g 2 contribute in the same way as explicit factors of m t . The algebraic reduction 
of two-loop tensor vacuum integrals proceeds along the same lines as described for the 
one-loop case above. The only difference is that tensors need not be totally symmetric 
beyond one loop. Finally, we are left with the scalar integral (R — S — 0), which can be 
easily calculated, 



r(m + n 2 + n 3 - n)T (f - m) r (m + n 2 - f ) r (m + n 3 - § ) 
r (f ) r(n 2 )r(n 3 )r(2n 1 + n 2 + n 3 - n) 



• (32) 



Now we identify the subgraph 7 with the top-quark loop. In this case the Taylor 
expansion in the integrand J 7 concerns the momenta that are external with respect to the 
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top-quark loop. Thus, all top-quark propagators P(q2 —p,m>t) are replaced by 



P(q 2 - p , mt ) = [( q2 - pf - m 2 ]' 1 = ]T(g 2 2 - m 2 )" 1 "^?^ - p 2 )', (33) 

where q2 is the loop momentum running through the top-quark loop. Note that p includes 
all external momenta of the process as well as the loop momentum q\ running through 
the internal gluon lines. The integration over q2 leads to one-loop vacuum integrals of the 
form fl28|), the calculation of which is described above. The integration over q\ involves 
only massless propagators to the first power. Since the q\ integration does not involve 
any m t -terms, it does not affect the power counting in m t at all. For one-gluon exchange 
the integral over qi is a massless tadpole, which vanishes in dimensional regularization. 
For two-gluon and three-gluon exchange the q± integration leads to the one-loop tensor 
integrals 

(27r^) 4 -™ r qi,^...qi, m 



B^...^ R (p) = — / d n q x 

m l J qiiqi+p) 

(2 tt//) 4 -" I ■ ■ ■ 

in 2 J 1 ql{qi+Pi) 2 {qi+P2) 2 ' 



c^M^) = ^/^ ^.T-;^-^ (34) 



The tensor integrals can again be recursively reduced to the corresponding scalar integrals 
in a fully algebraic manner. The relevant scalar integrals can be easily calculated and are 
given by 



B(p) 



4V \^ r( 2 -f)r(f-i) : 



pVo \-p 2 -iOJ T(n-2) 



2 



C (Pl,P2) 2 2 n 2 , = — ~2 ^ W ^ " 35 

pf=P2=°'P 2 =(pi-P2) 2 ^o p 2 y—p^ — iQj Y(n — 3) 

For p 2 = the integrals are zero in dimensional regularization. The case B(p) with p 2 = 
occurs, for instance, in graphs like Fig. §b. Therefore, we find that the contribution in 
the expansion (^) for which 7 is the top-quark loop is only non-vanishing in diagrams 
like Figs, ^c and d, where both external gluons are attached to the internal gluon lines. 

A few more "physical" remarks on the formally described algorithm for the asymp- 
totic expansion seem to be in order. The different contributions to a given Feynman graph 
that are associated with the subgraphs 7 in the expansion ([2B|) are directly related to the 
effective diagrams in the approach of the low-energy theorem described above. Shrinking 
the subgraph 7 to a point leads to the corresponding effective diagram, where the point 
arising from 7 is the point-like interaction of the effective Lagrangian. A non-vanishing 
contribution of a subgraph 7 7^ 7 loop requires that at least one external momentum pi 
passes through a massless propagator; otherwise the loop integral involving the mass- 
less propagators is zero. In other words, it is necessary that there exists a cut through 
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the diagram that passes exclusively massless lines. Therefore, only diagrams with such 
a "massless cut" can lead to contributions to effective diagrams in which an effective 
coupling appears in loops. But all diagrams in general contribute to tree-like effective 
diagrams, which result from shrinking the complete loop part 7i 00 p to a point. 

Finally, we mention that 75, which appears in the case of pseudoscalar Higgs bosons, 



is treated according to the prescription of 't Hooft and Veltman [25]. Technically we sub- 
stitute 75 by ie AtlAt2M3/14 7 Ml 7^ 2 7^3 7^/4! before the evaluation of the Dirac trace so that the 
actual trace calculation can be carried out for usual n-dimensional Dirac matrices. The 
correct projection of the trace result on the physical four-dimensional space is achieved 
upon the contraction with the four- dimensional e-tensor. In this approach, all loop in- 
tegrations can be carried out before the contraction with e, i.e. n-dimensional momenta 
can be used. Note, however, that the contraction with e necessarily occurs before the 
integration over the phase space of the radiated parton in the real corrections, i.e. one 
has to take care of the four-dimensionality of e there. Moreover, an additional spurious 



counter term has to be added to the Att vertex (see also |L4|, [18]]). 
3.4 Drell— Yan-like processes 




Figure 5: Diagram contributing to qq — > Ah, AH at lowest order. 



Pairs of scalar and pseudoscalar Higgs bosons can also be produced in qq collisons via 
s-channel Z-boson exchange, see Fig. [5[ At LO the partonic cross sections are given by 

a LO (qq - Ah, AH) = X ZAhjZAH ^-( Vq + 0ff ) _______ , (36) 

where Q 2 denotes the partonic cm. energy squared, and v q = 2I 3q — 4e q sin 2 Q w , a q = 2I 3q 
are the vectorial and axial charges of the initial-state quarks. The QCD corrections 
coincide analytically with the QCD corrections to the Drell- Yan process qq -—>■ Z* , if 
squarks and gluinos are heavy, so that their contributions can be neglected. Thus the 
NLO cross section can be expressed as: 

0"nloG°P A h, AH + X) = a LO + Aa qq + Aa qg , 
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c^lo = / dT^2——cr LO (Q 2 = Ts), 

Jtq a UT 



Act, 



71 



r 1 dC qq f 1 

/ dr^2 —7— / dz cr LO (Q 2 = rzs) LO qq -(z), 

Jtq „ CLT Jtq/t 



Act, 



a s (n) 



'/</ 



rl dC q9 f 1 

/ dT J2~] — / dz a LO (Q 2 = rzs) uj g g(z), (37) 



q,q 



with the coefficient functions |26 



-P M , 2 )lo g ^ + i{2[C 2 -2 Wl - 2 ) + 4(Mi^) 



-2(1 + ^) log(l-z) 



1 / M 2 

<»M = -2 P -W l0g ((T3^ 



)+^{l + 6,-7, 2 } 



(38) 



The Altarelli-Parisi splitting functions P q q and P qg are given by 




^K + (l-^) 2 }- 



(39) 



4 Results 

4.1 Standard Model 

The analysis for SM Higgs-boson pair production has been carried out for the LHC cm. 
energy y/s = 14 TeV. The top-quark and bottom-quark masses have been chosen as 
m t = 175 GeV and m b = 5 GeV. We have adopted the CTEQ4L and CTEQ4M |7 



parton densities for the LO and NLO cross sections, respectively, corresponding to the 
QCD parameters A5 = 181 MeV and A^ s = 202 MeV. Since our results have been 
obtained in the heavy-quark limit, they can be expected to be reliable for M# < 200 GeV, 
based on the experience from single-Higgs production via gluon fusion. 

In order to investigate the size of the QCD corrections, we define the K factor K = 
cnlo/clo as the ratio of the NLO and LO cross sections, where the parton densities and 
the strong coupling a s are taken at NLO and LO, respectively. The K factor for SM- 
Higgs pair production is presented in Fig. || as a function of the Higgs mass Mh and shows 
little variation with Mh- It ranges between about 1.9 and 2.0, thus enhancing the LO 
cross section significantly. Moreover, the broken lines of Fig. |6| show the contributions of 
the individual terms of Eq. (|20|) . It can be clearly inferred that analogously to the case of 
single-Higgs production the (infrared-regularized) virtual and real corrections originating 
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Figure 6: K factors of the QCD-corrected gluon-fusion SM cross section a(pp — > HH + 
X) at the LHC with cm. energy ^/s = 14 TeV. The dashed lines show the individual 



contributions of the four terms of the QCD corrections given in Eq. (\2ui). 
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Figure 7: The renormalization and factorization scale dependence of the SM-Higgs pair- 
production cross section at LO and NLO for a Higgs mass M H = 200 GeV. 
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from gg intial states dominate the QCD corrections, while the gq and qq initial states do 
not exceed about 5% in total. We note that the values shown for K vin , K gg , K gq , and 
K qq do not exactly add up to K tot — 1, since the individual contributions are obtained by 
taking NLO parton densities and the strong coupling in the NLO cross sections, but LO 
quantities in the LO cross sections consistently. 

In order to investigate the reliability of our results, the residual dependence of the 
cross section on the renormalization and factorization scales is shown in Fig. |7| for a 
Higgs mass of Mh = 200 GeV at LO and NLO. The scale dependence significantly 
decreases compared with the LO result by including the QCD corrections. However, the 
still monotonic decrease of the NLO cross section with increasing scales signals the need 
for further improvements. Nevertheless, the theoretical uncertainty of the NLO result can 
be estimated to about 20% from the residual scale dependence. 

The final results for the total Higgs pair-production cross sections at LO and NLO 
are presented in Fig. [8] as a function of the Higgs mass Mh- We recall that the full 
mass dependence of the LO form factors is included in all numerical evaluations, i.e. only 
the relative NLO correction is treated in the heavy-quark limit, in order to increase the 
validity of our results. For Mh ^ 200 GeV the cross section exceeds 10 fb, leading to 
more than about 3000 events at the LHC, once the anticipated integrated luminosity 
/ Cdt = 3 x 10 5 pb -1 is reached^. The typical signatures for Higgs-boson pair production 
are bbbb and bbr + r~ final states for Mh ^ 140 GeV. For Mh ^ 160 GeV, the Higgs-boson 
pairs decay predominantly into four vector bosons. 

Finally, we present the sensitivity of the SM-Higgs pair-production cross section to 
the trilinear Higgs coupling in Fig. |9], which shows the ratio of the total cross section 
with a non-standard trilinear coupling A and the SM one as a function of the trilinear 
coupling, varying in units of the SM coupling. The cross section becomes significantly 
larger for smaller trilinear couplings, so that this process may serve as a possibility to 
measure this coupling and test the SM Higgs sector, if the signal can be extracted from 
the QCD background. 

The trilinear Higgs coupling can potentially be measured at high-energy e + e~ colliders 



28|]. Since the rate for multiple Higgs-boson production at e + e~ colliders is small, high 
luminosities are required for a sufficient number of events in order to probe this coupling. 
However, the backgrounds will be well under control. 



4.2 Minimal supersymmetric extension 

Promising MSSM-Higgs pair-production processes at the LHC are hh, hH, HH, hA and 
HA production with sizeable ranges in the parameter space where the cross sections exceed 
10 fb [IT]]. The gg — > hh process can be used to cover a part of the MSSM parameter space 
for small tg/3 via its decay modes into 6677 final states ]29[ in the region where this process 
is dominated by resonant gg — > H — > hh production. Since our calculation is expected 
to be valid for small values of tg/3 (where the 6-quark contribution can be neglected) and 

3 The cross section at the Tevatron is less than 0.2 fb for y/s = 2 TeV. 
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Figure 8: SM prediction of the Higgs-boson pair-production cross section at LO and NLO 
as a function of the Higgs mass M H . 




Figure 9: Ratio of the Higgs-boson pair-production cross section at NLO with a non- 
standard trilinear coupling A and the SM cross section as a function of the trilinear cou- 
pling in units of the SM one. 
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Figure 10: K factors of the QCD- corrected gluon-fusion and Drell-Yan like cross sections 
a(pp — > 4>i4>2 + X) at the LHC with cm. energy y/s = 14 TeV. 



Higgs masses below the tt threshold, it may be assumed to reliably approximate the cross 
sections of light-scalar-Higgs pair-production in particular. Generally, in the MSSM, we 
expect our results to be valid for Ma ^ 200 GeV, while for larger pseudoscalar masses at 
least the heavy s-channel Higgs particles become too heavy with respect to the top-quark 
mass. 

In our analysis we have included the MSSM Higgs masses and couplings at the two- 
loop level in the effective-potential approach ||. Moreover, we have included all available 
higher-order corrections to the total MSSM Higgs decay widths [[JIJ, which appear in the 
s-channel Higgs propagators of the triangular loop contributions. 

The K factors for all Higgs-boson pair-production processes for tg/3 = 3 are presented 
in Fig. |10] as a function of the pseudoscalar Higgs mass M&. The numbers at the top are 
the corresponding values of the Higgs masses and Mh- The total K factors increase 
the Drell-Yan like cross sections by about 30% and the gg production cross sections 
by about 60-100%, thus signalling the importance of the QCD corrections. These are 
dominated by soft and collinear gluon radiation from gg and qq initial states, similar to 
the SM case discussed above in detail. The sharp decrease of the K factor for gg — > hh 
at Ma ~ 175 GeV originates from the resonance contribution gg —>■ H — > hh, which is 
kinematically forbidden below this mass and allowed above. The K factor of resonant 
single-Higgs production is smaller than the one of continuum hh production. 

In spite of the large size of the NLO contributions, the scale dependence significantly 
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decreases as can be inferred from Fig. 11 , which presents the scale dependence of the cross 
section a(pp — > hh + X) at LO and NLO for a light-scalar-Higgs mass = 95.5 GeV, 
corresponding to Ma = 200 GeV. Thus, the QCD-corrected results turn out to be reliable 
within about ±20%. The total NLO cross sections for the processes gg — > hh, hH, HH 
are presented in Fig. [12] as a function of the pseudoscalar Higgs mass Ma- The sharp 
increase of the hh cross section at Ma ~ 175 GeV is due to the fact that the resonant 
gg — > H —>■ hh process opens up above this mass value, while it is kinematically forbidden 
below. There are sizeable regions where the cross sections exceed a level of 10 fb. This 
is in contrast to AA production, the cross section of which is smaller than 1 fb in the 
considered parameter space. 

The total cross sections for the processes pp — > hA, HA and their individual contri- 
butions from the gg and qq initiated processes are presented in Fig. [13] as a function 
of the pseudoscalar mass Ma for tg/3 = 3. While for HA-p&ir production gluon fusion 
gg — > HA is always suppressed against the Drell-Yan like process qq — ► HA, both cor- 
responding processes are competitive for the light scalar Higgs particle h. Especially for 
smaller masses Ma the total cross sections for pp — > hA, HA exceed 10 fb and thus po- 
tentially provide the possibility to detect these processes. Note, however, that only the 
gluon-fusion processes are sensitive to the trilinear Higgs couplings. 

Finally, it should be noted, that in all neutral Higgs-pair processes the residual theo- 
retical uncertainties reduce to a level of 20% after including the QCD corrections. 

The corresponding NLO cross section for hh production at the Tevatron with cm. 
energy y/s = 2 TeV is presented in Fig. ||as a function of the pseudoscalar mass Ma- 
For Ma ^ 175 GeV it exceeds 10 fb, since the resonant gg —>■ H —>■ hh is kinematically 
possible and dominant in this mass region. The signatures of this process are bbbb and 
bbr + T~ final states. 



Trilinear MSSM Higgs couplings can be measured at high energy e + e~ colliders [pl|j , 
if sufficient numbers of signal events for Higgs pair production will be produced. At the 
necessary high luminosities background processes do not cause any problem. 



5 Conclusions 

We have presented a complete calculation of the two-loop QCD corrections to neutral- 
Higgs pair production at the LHC via gluon fusion in the limit of a heavy top quark. This 
approximation is at least reliable if the invariant mass of the produced Higgs-boson pair 
is below the tt threshold of the mediating top-quark loops. We have analyzed the results 
within the SM and the MSSM and found large increases of the cross sections by about 
60-100%. The QCD corrections to the associated production processes qq — > hA, HA 
coincide with those to the Drell-Yan process qq — > Z* , thus increasing the cross sections 
by about 30%. The QCD corrections stabilize the theoretical predictions compared with 
the LO results, which exhibit large theoretical uncertainties. After including the QCD 
corrections, the remaining theoretical uncertainties are reduced to a level of about 20%. 
Except for AA-pair production, all Higgs-pair production cross sections of the SM and 



21 



5x 1(T 



2x 10 



10 



I I I I 1 


111,1111,1111,111 
G(pp -> hh+X) [fb] 


1,11.1,1111,1111 

^ = M 2 = ^ 2 Q 2 




Vs=14 TeV 


m t = 175 GeV 






tg|3 = 3 


\ \. 
\ \ 
\ 

\ 




M h = 95.5 GeV 


\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 










LO 


, , , , 1 


, , , i , , , , i , , , , i , , , 


, i , , . , i , , , , i , , , , 



0.5 



1.5 



2.5 



3.5 



Figure 11: The renormalization and factorization scale dependence of the Higgs pair- 
production cross section a{pp — > hh + X) at LO and NLO for tgj3 = 3 and a Higgs mass 
M h = 95.5 GeV (M A = 200 GeV). 
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Figure 14: Scalar Higgs-boson pair-production cross section a(pp — > hh + X) at NLO 
as a function of the pseudoscalar Higgs mass Ma for tgf3 = 3 at the Tevatron pp collider 
with cm. energy y/s = 2 TeV. 

the MSSM exceed 10 fb at LHC energies in certain regions of the parameter spaces. As 
soon as the H hh channel opens, the cross section for pp hh + X reaches even 10 3 fb 
at the LHC and is still larger than 10 fb at the Tevatron for tg/3 = 3. 

Moreover, we have shown that the NLO prediction of Higgs-boson pair production is 
sensitive to a deviation of the trilinear Higgs coupling from its SM value, rendering this 
process well suited for studying this coupling in pp collisions. 
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